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ABSTRACT 

We present the results of three separate searches for Hi 21-cm absorption in a total of 
twelve damped Lyman-a absorption systems (DLAs) and sub-DLAs over the redshift range 
Znbs — 0.86 — 3.37. We find no absorption in the five systems for which we obtain reason- 
able sensitivities and add the results to those of other recent surveys in order to investigate 
factors which could have an effect on the detection rate: We provide evidence that the mix 
of spin temperature/covering factor ratios seen at low redshift may also exist at high redshift, 
with a correlation between the 2 1 -cm line strength and the total neutral hydrogen column den- 
sity, indicating a roughly constant spin temperature/covering factor ratio for all of the DLAs 
searched. Also, by considering the geometry of a flat expanding Universe together with the 
projected sizes of the background radio emission regions, we find, for the detections, that the 
21-cm line strength is correlated with the size of the absorber. For the non-detections it is 
apparent that larger absorbers (covering factors) are required in order to exhibit 21-cm ab- 
sorption, particularly if these DLAs do not arise in spiral galaxies. We also suggest that the 
recent z a b s = 2.3 detection towards TXS 031 1+430 arises in a spiral galaxy, but on the basis 
of a large absorption cross-section and high metallicity, rather than a low spin temperature 
(lYorketal.ll2007h . 

Key words: quasars: absorption lines - cosmology: observations - galaxies: high redshift - 
galaxies: ISM - radio lines: galaxies 



1 INTRODUCTION 

Damped Lyman-a absorption systems (DLAs) are believed to be 
the precursors of modern day galaxies, containi ng at least 80% 
of the neutral gas mass density of the Universe JProchaska et al.l 
2005). As the name suggests, DLAs are identified through their 
heavily damped absorption features, due to the large columns of 
neutral hydrogen (Nui ^ 2 x 10 20 cm -2 ) through which the back- 
ground quasar is viewed. The 21-cm spin-flip transition is of inter- 
est since it traces the cool component of the gas, with the compar- 
ison of the 21-cm absorption strength to the total neutral hydrogen 
column giving the gas spin temperature (T sp i n ) for a fully absorbed 
emission region (/ = 1) [see Equ.[T] Sect. 13.11 . Searches for 21-cm 
absorption in DLAs exhibit a ~ 50% detection rate, the detections 
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occuring predominately at low redshift (z a bs ~ 1), suggesting an 
increase in the T sp i n /f ratio with redshift, which could be due to 
higher spin temperatures and/or lower covering factors. 

In order to shed light on which is the predominant factor, we 
have undertaken searches with the Green Bank (GBT) and Giant 
Metrewave Radio (GMRT) Telescopes. In particular, we aim to: 

(i) Test the hypothesis o flCurranetai]j2005h that 21-cm absorp- 
tion should be readily detectable at high redshift, despite the large 
Tspin// ratios, by searching for 21-cm absorption in DLAs which 
lie towards compact radio sources. 

(ii) Test the hypothesis of ICurran eT al. (2007b) that the spin 
temperature does not continue to increase with redshift, by search- 
ing in previously unsearched high redshift (z a b s ~ 3.2) DLAs. 

(iii) Test the line stren gth-metallici t y corre lation (T sp m//- 
[M/H] anti-correlation) of Cur ran et alj 112007b). which suggests 
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that several, currently undetected, DLAs should be readily de- 
tectable in 21 -cm absorption. 

We have observed a dozen DLAs and sub-DLAs and in this paper 
we present and discuss the results in the context of the above issues. 



2 SAMPLE SELECTION AND OBSERVATIONS 
2.1 GBT 2006 observations 

2.1.1 Motivation 



Prior to the recent hi gh redshift detections ("KanekaretaD[2006, 
120071 : lYork et alj 120071 discussed below), all of the DLAs de- 
tected in 21 -cm absorption were at redshifts of z a b s ~ 2, where 
there is also an approximate ly equal number of non-detections. 
iKanekar & Chengalun J2003h . interpreted this as the low redshift 
DLAs having a mix of spin temperatures, whereas the high redshift 
absorbers have exclusively high spin temperatures. However, all of 
the Zabs ~ 1 absorbers have large DLA-to-QSO angular diameter 
distance ratios (-DAdla/-D^4qso w 1), whereas the z^hs ~ 1 ab- 
sorbers have a mix of ratios, due to the ge ometry effects of a flat 
expanding Universe dCurran & W ebb 2006): Those with low ratios 
are almost always detected in 21 -cm absorption, whereas those with 
high ratios tend not t o be detected, which is how th e "spin temper- 
ature" distribution of Kanekar & Chengalurj d2003l) is segregated. 
This strongly suggests, on the grounds of geometry alone and not 
allowing for various absorber and emitter sizes, that the covering 
factor, rather than the spin temperature, is the crucial criterion in 
determining the detectability of 21 -cm absorption, where an ab- 
sorber at high redshift will always cover the background quasar 
much less effectively than an equivalent absorbe r-quasar system 
at low redshift. Furthermore, ICurran et alj J2005h had previously 
suggested that the non-detections at high redshift could be due to 
low covering factors, as opposed to elevated spin temperatures, and 
so 21 -cm absorption should be readily detectable in high redshift 
DLAs towards very compact radio sources. 

With the Green Bank Telescope we therefore embarked upon 
a survey for Hi 21 -cm absorption in the known radio-illuminated 
DLAs yet to be searched, prioritised by a small radio source size. 
For this we used the Sloan Digital Sky Survey Dampe d Lya Sur- 
vey Data Release 1 dProchaska & He rbert-Fort 2004) combined 
with the previously unsearched systems dCurran et alj|2002h . giv- 
ing 57 DLAs occulting radio-loud (S iS 0.1 Jy) quasars. We then 
short-listed the DLAs and sub-DLAs in which the 21 -cm transition 
was redshifted into the Prime Focus 342 and 450 MHz bands and 
which were located towards the four most compact radio sources - 
B3 0758+475 (SDSS J080137.67+472528.2), PKS 1402+044, PKS 
B1418-064 and PKS 2136+141 (Tabled where we also give the 
radio source sizes for the other objects searched). 



Table 1. The radio source sizes at ug of the QSOs illuminating the DLAs 
searched (Table [2) used for selecting the GBT 2006 targets. These were 
obtained from the highest resolution radio images available, closest in fre- 
quency to the redshifted 21-cm values (yg), although 0454+ 039, 0528-250 
and 1418-064 have since been imaged at closer frequencies (Kanekar et al. 



2009a). For the remainder of the DLA sample refer to Curran et al. (2005 
2007bi) . 



QSO 


e QSO [arc-sec] 


Ug [GHz] 


Reference 


0454+039 


< 0.4 


2.3 


Beaslev et al. (2002) 


0528-250 


0.010 


5.0 


Fomalont et al. (2000) 


0758+475 


0.91 


1.4 


FIRST 


J08 16+4823 


0.35 


1.4 


FIRST 


0957+561A 


11.36 


5.0 


Harvanek et al. (1997) 


1402+044 


1.35 


1.4 


FIRST 


1418-064 


1.27 


1.4 


FIRST 


2136+141 


0.002 


5.0 


Gurvits et al. (1999) 



FIRST: The Very Large Array's Faint Images of the Radio Sky at Twenty 
Centimetres. 
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Figure 1. GBT 336 MHz band spectrum towards 0758+475 shown at a ve- 
locity resolution of 10 km s _1 . The arrow in this and the following spectra 
(Figs. l2"lto l 101 indicates the expected frequency of the absorption. 



spacing of the GBT spectrometer gave a resolution of 0.34 km s~ 1 . 
The system temperature was typically T sys ~ 120 K and RFI was 
not too severe, allowing 6.0 hours of data to be retained. Although 
most spectral baselines were flat with the odd intermittent spike, the 
measured fluxes were found to vary considerably between scans 
(from ~ — 1 to ~ 1 Jy), not permitting us to determine the flux 
density of this source (Fig.[TJl. 

1402+044 was observed at a central frequency of 385 MHz 
on the 30th of November and 1st of December 2006, for a total of 
1.7 hours. As per previous WSRT observations towards this object 
(described in lCurran et alj|2007bl) . the relatively wide band of the 
GBT spectrometer could cover the z = 2.688 (38 5.14 MHz) and 
_2.713 (382.55 MHz) candidate sub-DLAs of iTurnshek et al.1 
, as well as the z = 2.708 DLA (383.07 MHz). At 385 MHz, 
the GBT spectrometer gives a spectral resolution of 0.30 km s~ x . 



Typically, the system temperatures were T SJ 



110 K and the 



removal of RFI dominated scans left a total observing time of 1 .4 
hours (Fig.[2j. 



2.1.2 Observations and data reduction 

The observations were performed with the GBT PF1 342 and 
450 MHz receivers in position switching mode (level-3 sampling, 
one spectral window and dual polarisation). The backend was either 
the GBT spectrometer using 32 768 lags over a 12.5 MHz band, 
giving a channel spacing of 381.47 Hz, or when RFI appeared to 
be particularly bad, the spectral processor, using 1024 lags over a 
10 MHz band, giving a channel spacing of 9.766 kHz. 

0758+475 was observed on the 1st of December 2006 for a 
total of 7 hours. At a central frequency of 336 MHz, the channel 
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Figure 2. GBT 385 MHz band spectrum towards 1402+044 shown at a 
velocity resolution of 10 km s~ x . 
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Figure 3. GBT 408 MHz band spectrum towards 1402+044 harming 
smoothed to a velocity resolution of 14 km s . 
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Figure 6. GMRT 336 MHz band spectrum towards J0801+472 (0758+475). 
The velocity resolution is 3.5 km s —1 . 
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Figure 4. GBT 320 MHz band spectrum towards 1418-064 hanning 
smoothed to a velocity resolution of 18 km s — 1 . 



The z = 2.485 absorber was observed using the spectral pro- 
cessor, tuned to a central frequency of 408 MHz, on the 6th and 8th 
of January 2007 for a total of 7.2 hours, of which 6. 1 hours were 
retained (T sys ~ 70 K). The channel width was 9.766 kHz, giving 
a spectral resolution of 7.2 km s _1 at 408 MHz (Fig. [3}. 

1418-064 was observed using the spectral processor, tuned to 
a central frequency of 320 MHz, on the 22nd of March 2007 for 
a total of 2.4 hours, of which 1.9 were retained (T sys ~ 80 K). 
As per 0758+475, although the bandpasses were flat, the measured 
flux density varied considerably between scans not permitting us 
to determine a value for this (Fig. [4). The spectral resolution was 
9.2 kms" 1 at 319 MHz. 

2136+141 was observed using the spectral processor, tuned 
to a central frequency of 454 MHz, on the 5th and 7th of January 
2007. Originally, 12 hours were planned to search for 21 -cm ab- 
sorption in the two (admittedly weak, Tabled absorbers towards 
this compact quasar (Table [TJ, but due to severe RFI over most of 
this band, we only observed for 2.0 hours (where T sya ~ 65 K). At 
redshifts of z = 2.118 and 2.134, the 10 MHz band centered on 
454 MHz covered both 455.6 and 453.2 MHz at a resolution of 6.5 
km s~ , with the RFI dominating 451 to 455 MHz. This allowed us 
to only obtain a limit for the lower redshift sub-DLA, which itself 
may not be reliable (Fig. [5}. 

All of the data were reduced using the GBTIDL package. Af- 
ter flagging, both polarisations were combined, linear (1st order) 
baselines removed and the resulting data Gaussian smoothed, via 
the gsmooth task, to a resolution of 3 km s - . The results are 
summarised in Table|2] 
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Figure 5. GBT 454 MHz band spectrum towards 2136+141 hanning 
smoothed to a velocity resolution of 13 km s . 



2.2 GMRT 2008 observations 

2.2.1 Motivation 

For the GMRT sample the intent was to find new high redshift 
"end points" with which to verify the flattening off of T ap i n // at 
« 2000 K at z abs £ 2 jCurran et al]|2007bh . To this end, we 
trawled the Sloa n Digital Sky Survey Damped Lya Survey Data 
Releases 1-3 JProchaska et al.l 120051) for DLAs occulting radio- 
loud quasars yet to be searched for 21 -cm absorption. This yielded 
two DLAs occulting quasars with flux densities in excess of 50 mjy 
at 1.4 GHz, J080137.6+472528.2 and J081618.9+482328.4. 



2.2.2 Observations and data reduction 

For all of the observations, all 30 antennae were used with the 325 
MHz receiver, where for an antenna gain of 0.32 Ky Jy _1 we ex- 
pected to reach an r.m.s. of o w 3 mjy per channel in each of the 
two polarisations. For the backend we used a correlator bandwidth 
of 0.5 MHz, giving a velocity resolution of « 4 km s" 1 and a cov- 
erage of « ±230 km s" 1 (Az « ±0.003). 

J0801+472 (which is 0758+475, Sectl2~H was observed at a 
central frequency of 336.366 MHz on the 7th and 8th January 2008. 
3C 147 was used for the bandpass calibration and 0834+555 was 
observed every 0.6 hours to correct for delays. In general, RFI was 
minimal and the data were excellent with all 30 antennae operating, 
although pairs involving antennae 13 & 28 were flagged due to poor 
phase stability. On the first night, the target was observed for a total 
of 7.54 hours, of which only 0.14 hours required flagging and 350 
good baseline pairs were retained. On the second night, the target 
was observed for a total of 7.78 hours of which 7.64 hours were 
kept. Again the data were of high quality, giving 400 good baseline 
pairs. After analysing each polarisation of each night separately, a 
combined image was produced using the invert task in MIRIAD. 
The source was unresolved by the 13" x 11" beam. 

J0816+4823 was observed at a central frequency of 320.214 
MHz for a total of 7.8 hours on 16th March 2008. 3C48 was 
used for the bandpass calibration and again the phase calibrator 
0834+555 was observed every 0.6 hours during the observation. 
After flagging the non-operating antennae, 8 & 28, as well as an- 
tenna 6, on which the phase stability was poor, RFI appeared mini- 
mal and only a little time dependent interference had to be removed 
(leaving 7.65 hours of data). However, perhaps due to its low flux 
density, an accurate image of the quasar could not be produced. 
We therefore flagged noisy antenna pairs until the averaged visi- 
bilities gave an r.m.s. noise level close to the « 10 mJy expected. 
Averaging the « 250 remaining baselines in each of the polarisa- 
tions, gave a similar flux density to extracting this from the cube, 
i.e. 85 mjy at 0.32 MHz, which compares well with the 78 mjy 
at 1.4 GHz from the The Very Large Array's Faint Images of the 
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Figure 7. GMRT 320 MHz band spectrum towards J0816+4823. The ve- Figure 9. GBT spectrum towards 0528-250 over the relatively RFI free 

locity resolution is 3.7 km s — . 372.53-373.20 MHz, shown at a velocity resolution of 3 km s . 
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Figure 8. GBT spectrum towards 0454+039, shown at a velocity resolution 
of 10 km s _1 . 
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Figure 10. GBT spectrum towards 0957+561A over the relatively RFI free 
590.6-595.4 MHz, shown at a velocity resolution of 10 km s . 



Radio Sky at Twenty Centimetres (FIRST). Again, the results are 
summarised in Table[2] 

2.3 GBT 2008 observations 

2.3.1 Motivation 

According to the 21 -cm line s trength-metallicity correlation 
(T S pi n // - [M/H] anti-correlation. ICurran et al]|2007bh . there are 
three previously searched DLAs (PKS 0454+039, PKS 0528-250 
& [HB89] 0957+561) for which the spin temperature/covering fac- 
tor ratio should be T sp i n / / < 1000 K, with the previous limits 
possibly being close to those required to detect 21-cm absorption. 
We therefore applied for Green Bank Telescope time in order to 
improve these limits to the point where 21-cm absorption would 
be detected, although a significant improvement was only managed 
for 0528-250 (see Sects.l3~Tl&l4~2T 

2.3.2 Obsen : ations and data reduction 

As per the 2006 run, we used the GBT Prime Focus Receivers in po- 
sition switching mode, with the GBT spectrometer split into 32 768 
lags over a 12.5 MHz band, giving a channel spacing of 763 Hz, in 
dual polarisation. 

0454+039 was observed for a total of 2.9 hours over several 
runs from 28th August to 13th September 2008, at a central fre- 
quency of 763.8 MHz, giving a spectral resolution of 0.3 km s _1 . 
System temperatures ranged from T sys ~ 30 to 40 K over these 
dates and the overall RFI was minimal, except over the frequency 
range of interest, 763.6 to 763.8 MHz, where a blooming great 
spike appeared in most of the spectra. This led to significant flag- 
ging, leaving 1.2 hours of good data, of which the average spectrum 
is shown in Fig. [8] 

0528-250 was observed 2nd and 18th February 2008 for a 
total of 7.4 hours. At a central frequency of 327.7 MHz, giving 
a spectral resolution of 0.7 km s - . The system temperature was 
typically T sys ~ 60 K and RFI was only severe over 0.64 hours of 
the observations. However, besides this there were intermittent RFI 
spikes over the whole band, which only permit r.m.s. estimates over 
the intervening clean regions. Therefore the limit quoted in Table|2] 



is for 372.53 to 373.20 MHz (z = 2.8060 - 2.8116), close to that 
expected from the redshift of the DLA (Fig. [9). 

0957+561A was observed for a total of 3.6 hours over several 
runs from 7th August to 8th November 2008, at a central frequency 
of 594 MHz, giving a spectral resolution of 0.4 kms^ 1 . System 
temperatures ranged from T sys ~ 30 to 60 K over these dates and 
RFI was severe. This led to extensive flagging of bad data, leav- 
ing only 14 minutes of "good data", of which the average spec- 
trum is shown in Fig. \10\ Although this flagging removed all of 
the spikes close to 594 MHz, a strong ringing was still present in 
the unsmoothed data and, as seen from the spectrum, no reasonable 
limit can be obtained from these observation^. 



3 RESULTS 

3.1 Observational results 

Typically the illumination of the optical (Lyman-a) and radio (21- 
cm) absorbers arise from quite different sources sizes, these be- 
ing < 1 pc for the former, but in excess of 100 pc for the latte^ 
(Wolfe et al. 2005 and references therein). Thus, it is possible that 
the optical and radio observations coul d trace different sig ht-lines 
(as is probably the case for 1622+238. ICurran et al.f 2007a). How- 
ever, if the 21-cm and Lyman-a absorption do arise in the same 
cloud complexes dPickev & Lockman| [l990), the column density 
[cm -2 ] of the absorbing gas in a homogeneous cloud is related to 
the velocity integrate d optical depth, wh ere r = — In (l — -jg), 
of the 2 1 -cm line via jWolfe & Burbidgdl 19751) : 

Nm = 1.823 x 10 18 T spin J t dv , (1) 

where T sp m [K] is the mean harmonic spin temperature of the 
gas, a is the depth of the line (or r.m.s. noise in the case of a 
non-detection) and S & f the flux density and covering factor of 
the background continuum source, respectively. In the optically 

1 Note that the unflagged data give a flux density of f» 1 Jy, although the 
spectrum is plagued by RFI spikes close to the tar get frequency 

Even over 10 kpc in some cases (see figure 4 o flCurran et alj2005h . 
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Table 2. Our search results. z abs is the redshift of the DLA giving the observed 21-cm frequency, f b s - R UN gives the observing run (Sects. |2~T1 
to 12.31 . S is the flux density at z^ b s , o" is the r.m.s. noise level reached per Av channel [km s ], where these are > 3 km s -1 ; for the GBT 
observations the resolutions are typically much finer (Sect. |2~TV Since the measured r.m.s. noise is dependent upon the spectral resolution, as usual 
jCurran et al . 2005, 2007b), for the non-detections the 3<r upper limits of T pc ^ at a velocity resolution of 3 km s — 1 are quoted. The total neutral 
hydrogen column density, TVjji [cm -2 ], is given with the corresponding reference, which yields the quoted spin temperature/covering factor ratio, 
via the estimated FWHM [km s — +] for the non-detections (see main text). 



QSO 


Zabs 


u ohs [MHz] 


RUN 


S[Jy] 


a [mJy] 


Ac 


^peak 


logic) N m 


Ref. 


FWHM 


T spin // [K] 


0454+039 


0.8596 


763.8 


GBT08 


0.70 


< 4.6 


3.0 


< 0.020 


20.7 


S95 


19 


> 690 


0528-250 


2.811 


372.7 


GBT08 


0.561 


< 6.9 a 


3.0 


< 0.037 


21.3 


L96 


40 


> 700 


J080 144725* 


3.2228 


336.4 


GBT06 




< 6.1 


3.0 




20.7 


P05 












GMRT 


0.24 


< 2.9 


3.5 


< 0.039 






20+ 


> 330 


J08 16+4823 


3.4358 


320.2 


GMRT 


0.08 


< 6.3 


3.7 


< 0.29 


20.8 


P05 


20+ 


> 60 


0957+561A 


1.391 


594.0 


GBT08 










20.3 


T93 


30 




1402+044 


2.713 


382.5 


GBT06 


1.41 


< 16 


3.0 


< 0.034 




T89 








2.708 


383.1 


GBT06 


1.41 


< 16 


3.0 


< 0.034 6 


20.9 


P04 


20+ 


> 610 




2.688 


385.1 


GBT06 


1.37 


< 20 c 


3.0 


< 0.044 




T89 








2.485 


407.6 


GBT06 


1.51 


< 4.8 


7.2 


< 0.015 


20.2 


T89 


20+ 


> 270 


1418-064 


3.449 


319.3 


GBT06 


0.374 d 


< 8.1 


9.2 


< 0.12 


20.4 


E01 


25 


> 40 


2136+141 


2.134 


453.2 


GBT06 










19.8 


T89 








2.118 


455.6 


GBT06 


0.320 


< 7.7 


6.4 


< 0.11 


19.8 


T89 


20+ 


> 10 



Notes: *0758+475 - observed with both the GBT and GMRT (Sects.[2j]&[272). +No Mg II equivalent width or metallicities found, s o 20 km a" 1 
assume d. "Over 372.53 to 373.20 MHz (z = 2.8060 - 2.8132). ^Double the sensitivity reached with the WSRT observations JCurran et al] 
l2007bl) . so we adopt this limit. c Subject to an RFI spike from 385.17 to 385.23 MHz (z = 2.6872 - 2.68 77). d Flux density could not be 
determined from our observations and so have adopted value for TXS 1418-064 at 0.374 Jy iDouglas et al Jl996l) . although this is 4'6" from PKS 
B 1 4 1 8-064 (for w hich S2.7GHZ = 0-38 Jy & S s.oOHz = 0-34 Jy) i t is wi thin the Texas survey beam . 

Refere nces: T89 -iTurnshek et all ll989|) T93 - ITurnshek & BohlinN 19931), S 95 - ISteidel et all jl995h , L96 - iLu et all ll996l) . E01 - lEilison et al] 
j200ll) . P04 - IProchaska & Herbert-FortM2004 . P05 - IProchaska et all <200"51) . 

of the absorbers and the limits for a further three are so poor as 
to be of little use. Of the remaining five, for 04 54+039 we have 
match ed the previous limit of r B p in // £ 600 K (Briggs & Wolfe 
[l983h . while for 0528-250 we have si gnificantly improved upon 
the T spin // > 100 K limit of lCarilli et all (l9960 This is the first 
search for 21-cm absorption at z a b s = 3.223 towards J080 1+4725, 
although, due to the relatively low flux density, the limit is only 
T spin // > 300 K0For the z abs = 2. 708 DLA towards 1 402+044, 
we have improved the previous limit dCurran et alj2007bl) by a fac- 
tor of two and for the z a b s = 2.485 sub-DLA, publish the first 
limit, although this is weak due to the low neutral hydrogen col- 
umn density. 



3.2 Update of the previous results 

Adding the new results and applying the estimated FWHMs, in 
Fig.QTJwe update the s pin temperatu r e/cove ring factor-absorption 
redshift distribution of ICurran et alj J2005h . From the histogram 
it is seen that the low redshift absorbers have a similar T S pin// 
distribution to those at high redshift (z a b s ^ 1-6) and in Table [3] 
we present the results of the statistical tests performed on the low 
and high redshift samples. These were done for redshift cuts of 
z CUT = 1.0, 1.6 ( the turnover in angular diame ter distance, Sect. 
|4TB and 2.0 (as Kanekar & Chengalur 2003). From these, we 



4 Based upon our 3 km s 1 resolution and 40 km s 1 line- width (Ta- 
ble|2). 

5 For this and the other GMRT target, J0816+4823, it is possible that 
the redshift coverage of z a b s ± 0.003 is not sufficient to cover uncer- 
tainties iritheabsorber redshift, which are determined from the Lyman-a 
line jProchaska et al. 2005). However, J0801+4725 has been searched over 
•Sabs = 3.30 ± 0.08 with the GBT (Fig. Q} and the J0816+4823 T spin /f 
limit is too poor in any case to contribute to the analysis. 



thin regime (<r <§C f-S), Equ. [TJ reduces to Nui = 1.823 x 
10 18 T *y" J ^ dv, which we apply to the entries in Table[2] where 
we summarise the results for the three sets of observations. 

For our and the previous non-detections we have no knowl- 
edge of what the FWHM of the 21-cm profile should be, thus 
hampe ring the determination of a lower limit for T sp m//. Previ- 
ously dCurran et al]|2005l,l2007bl) . we have used the » 20 km s -1 
mean of the detections to derive velocity integrated optical depth 
limits for the non-detections. An alternative to this is to esti- 
mate the FWHMs from the spin temperature, via the kinetic tem- 
perature according to T spi n w T kin < 22 x FWHM 2 (e.g. 
iKanekar & Chengalur! l2003h . although in the absence of detailed 
knowledge of the covering factor, it is impossible to assign a spe- 
cific spin temperature to any of the DLAs. In light of the findings of 
ICurran et 

ID d2007bj) , however, we can make educated guesses to 
the profile widths according to the FWHM- Mg II rest-frame equiv- 
alent width correlation for the detections. This gives FWHM~ 
13W r A 2796 (figure 6 of ICurran et afll2007bl) . and, where the Mgll 
equivalent widths are not available (generally at z a b s ~ 2.2), we 
estimate thes e from the metallicity via W^ 2 796 ~ 2.0 [M/H] +4 .0 
(figure 7 of ICurran et all l2007bL see also iMurphv et alj 12 007 fl 
Note that, although lGupta et aT. (2009): lKanekar et alj (2009b) find 
that the FWHM-W^ 2796 correlation does not generally hold for 
Mg II absorption systems, when only those which are likely to be 
D LAs (i.e. with an Mg I 2852 A equivalent width of W^ 2852 > 0.5 
A. lRao et alj|2006l) are added to the s ample, the c orrelation is still 
found to hold for all suspected DLAs ICurran 2009). 

Commenting on our observational results, due to the effects 
of RFI and weak radio fluxes, we cannot assign limits to four 

Applying the above conversions gives a range of 5 - 42 km s _1 for the 
non-detections, thus not drastically changing the values of log 1 g(T sp ; n //), 
shown in the following figures, from those using a FWHM of 20 km s~ 1 . 
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Figure 11. The spin temperature/covering factor ratio versus the absorption redshift for the DLAs searched in 21-cm absorption. The shapes represent the 
type of galaxy with which the DLA is associated: circle-unknown type, star-spiral, square-dwarf, triangle-LSB. The filled symbols represent the 21-cm 
detections and the unfilled symbols+arrows show the lower limits for the non-detections (all of these bar one, 0454+039 at z abs = 0.8596, have unknown 
host identifications ). The grey symbols s how the results since ICurran e t al. (2005) and the coloured symbols show the new results presented here [blue] 
as well as those of Kanekar et al. (2009b) [magenta]. The solid grey histogram shows the spin temperature/covering factor distribution for the low redshift 
( z abs < 1-6) absorbers and the outlined histogram the distribution for the high redshift (z a b s ^ 1-6) absorbers, where the limits to T sp i n // are treated as 
detections. 



Table 3. The statistics for and between the low and high redshift samples for various redshift cuts, where P gives the probability 
that the low and high redshift DLAs are drawn from the same sample and S gives the corresponding significan ce (assuming Gau ssian 
statistics). The values in the first four rows incorporate the limits to the non-detection s via the ASURV pa ckage llsobe et al. 1986) and 
for the Kolmogorov-Smirnov test the limits are treated as detections (as per Kanekar & Chengalur 2003). 

Z CUT = Z CUT = 2 CUT ~ 

^abs ^ 2 CUT ^abs ^ ^CUT ^abs ^ Z C\JT ^abs ^ ^CUT ^abs ^ ^CUT ^abs ^ 2 CUT 

Sample size, n 19 20 23 18 25 16 

Mean, log 10 (T apin //) 3.08 ±0.18 3.38 ±0.18 3.16 ±0.17 3.26 ±0.16 3.11 ±0.15 3.30 ± 0.20 

^Log-rank 0.0718 0.265 0.203 

■^Log-rank l-80o- 1.39ct 1.27a 

Pks 0.443 0.448 0.358 

Sks 0.767cr 0.759o- 0.919<x 



find slightly lower mean values of T sp i n // for the low redshift ab- 
sorbers, with log-rank and Kolmogorov-Smirnov tests indicating 
that there is little significance in the differences between the low 
and high redshift samples. Bear in mind, however, that the FWHMs 
for the non-detections have been estimated (Sect. 13. It , with some 
of these relying on a secondary correlation (FWHM oc W^ 2796 ~ 
2.0 [M/H] + 4.0), when the Mgll equivalent width is not avail- 
ably. Furthermore, the Kolmogorov-Smirnov test treats the upper 
limits as actual values and higher sensitivity limits at high redshift 
are required in order to confirm the similarities between the sam- 
ples apparent in Table [3] 

Performing a Kolmogorov-Smirnov test on the redshift dis- 
tributions between the detections and non-detections, thus requir- 
ing no assumptions, gives an 11% probability (significant at 1.6a) 



e Setting the FWHMs obtained from the metallicity to 20 km s _1 , how- 
ever, alters only the z a t, s z cut averages in Table [3] and these only 
slightly. Overall, there is very little change with 5L og — rank = 1-95 1 1-45 
& 1.45cr and S K S = 1-15 , 1.22 & 0.42cr, respectively. 



that these are drawn from the same sample. Additionally, the non- 
detections have a higher mean redshift - = 1.77 (a = 1.08), 
compared to z a bs = 1-13 (a = 0.96) for the detections. There- 
fore, we cannot rule out that the T sp i n // distribution differs be- 
tween the two redshift regimes. Although the relative contributions 
of T S pin and / are unknown, the fact that the largest differences 
are seen for z CUT = 1.0 (where [T sp i n //] 2 <i = 1200 K cf. 
[r ap in//] z^i = 2400 K) is consi stent with geometry ef fects giv- 
ing larger values of / at z a b s < 1 dCurran & Web bll2006| and Sect. 
[4X3. 

While we are discussing statistical trends over the whole sam- 
ple, it should be emphasised that the 21-cm absorption searches 
yield a very heterogeneous sample, with each survey having dif- 
ferent selection criteria (including the three here), resulting in a 
mix of detections and non-detections across various absorber mor- 
phologies (where these are known, Fig. lilt . However, presuming 
that the optical and radio absorption clouds are coincident, Equ. 
Q] accounts for differences in the neutral hydrogen column densi- 
ties and the flux densities o f the background sources (discussed in 
detail in ICurran et al]|2005t) . leaving only the T sp i n // degeneracy. 
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Figure 12. The velocity integrated optical depth of the 21 -cm absorption 
versus the total neutral hydrogen column density. The symbols are as per 
Fig.Gj] 



Naturally, various properties of the absorber and emitter, as well as 
the distance between them, could influence this ratio, issues which 
we discuss here. However, whether it be due to similar values of 
Tspin and f for all of the absorbers, or just the ratio being main- 
tained, there may be no large differences between the low and high 
redshift samples. Furthermore, since / ^ 1, the ordinate in Fig. II II 
represents the maximum possible spin temperatures for the detec- 
tions of which there are now several at z a b s ~ 2. This indicates 
that high redshift s do not necessarily entail hig h spin temperatures, 
as advocated by IKanekar & Chengaluj 1 20031) [although a larger 
fraction of warm gas jKanekar et al.ll2009ah does not exclude the 
possibility of cold clumps]. Case in points against exclusively high 
spin temperatures at high redshift are, 



(i) The 21 -cm detection at z abs = 2.289 towards 0311+430 
dYork et al]2007» , which has one of the lowest "spin temperatures" 
yet found, i.e. T spln /(/ = 1) » 140 K. 

(ii) The excitation temperature of T cx ~ 50 K derived from 
the H 2 transitions in the z abs = 2 626 DLA towards 08 12+32B, 
indicating the presence of cold gas Jjorgenson et~al]|2009h . 

(iii) The other 17 D LAs in which Eb has been detected 
dNoterdaeme et alj|2008l and references therein). These are gener- 
ally at z a b s ~ 1.7, although they may represe nt a specific class of 
absorber within th e general DLA population dCurran et alj|2004l ; 
Murp hy et al]|2004h . 



Lastly, ICurran et al. I d2005h also noted a correlation between 
the integrated 21 -cm optical depth and the Hi column density. 
Updating this for the DLAs (Fig. IT2T > and incorporating the upper 
limits (as above), there is a probability of P(t) = 0.0021 that the 
correlation occurs by chance. This corresponds to a 3.08u signif- 
icance, and since the gradient is oc //T sp i n , suggests that, on the 
whole, the ratio of spin temperature to covering factor is a constant 
for DLAs, although there is some considerable scatter. 



4 DISCUSSION 
4.1 Covering factors 

4.1.1 Angular diameter distances and other effects 

As mentioned in Sect. |2.1| Curran et alJ l feOOsT) suggested that, con- 
trary to the high redsh i ft equ als high spin temperature argument of 
IKanekar & Chengalud l2003h , 21 -cm absorption should be readily 
detectable at high redshift towards compact radio sources, where 
the covering facto r is maximised. At the t ime of writing up the re- 
sults presented in ICurran & Webbl (2006), this was confirmed by 
the detection of 21 -cm absorption in the hig hest redshift exampl e 
to date, at z abs = 2.347 towards 0438-436 jKanekar et alJl2006h . 
which has a backg round radio source size of only 0.039" at 5 GHz 
<TingavetalJll998h . Following this, since submitting the proposal 
for the 2006 GBT observations, 21-cm absorption has been con- 
firmed at yet higher redshift, z a b s = 3.386 towards 0201+113 
{IKanekar et alj|2007l see also Ide Bruvn et al.l fl996l : iBriggs et al 
I l997h . as well as at z abs = 2.289 towards 0311+430 jYork et al 
120071) . Again, both of these occult compa ct radio sources of < 

1990£|and 1.36" x 0.63", 



0.007" at 1.4 GHz (Stanghellin i et al 



respectively. The latter was noted by lYork et alj d2007l) to be un- 
usual in its low spin temperature of T sp m ^ 140 K at such a high 
redshift, although t his is not unusual in re gard to the finding of 
ICurran et alj fc005» : ICuiTan & Webbl j2006h . 

From Table Q] we see that five of our targets are towards 
radio source sizes of ^qso ~ 1" (re cently found to be .07" 
for 0454+039 and 0.17" for 0528-250, IKanekar et ai]|2009al) . al- 
though we have only meaningful T sp - m / f limits for three of these 
(0454+039, 0528-250 & 0758+475 - Table|2j. However, all of the 
GBT 2006 targe t s (Sec t. l2~Tt were selected before the arguments of 
Curran & Webb ( 2006) were formulated, and so all of the 2006 tar- 
gets (and two of the three just mentioned), are all at high redshift, 
and will therefore be disadvantaged with respect to the DLA-to- 
QSO angular diameter distance ratios. 

Showing these in Fig. [T3] (top), we see that all but one of 
our target DLAs (for which we have limits), are indeed at angu- 
lar diameter distance ratios of DAola/D Aqso ~ 1, due to their 
high redshift selection. For a given absorber and radio source size, 
this will disadvantage the effective covering factor (possibly giving 

, Sect. l3.2b , although, 



the observed [T sp i B /f]^i w 2 [T spin /f] z<1 
as mentioned above, the three new high redshift detections occult 
compact radio sources, as shown in the bottom panel of F ig. 1131 
where we use the source sizes compiled i n | Curran et"aD d2005t) up- 
dated with those given in IKanekar et ail d2009iH 8 l . However, all of 
our non-detections (which are predominately at z a b s ~ 2.4), also 
occult similar radio source sizes. Of these, three of the sizes are 
not as reliable (orange symbols where 3 < Vsize/Vabs < 10, cf. 
green where fsme/vabs ^ 3), although as seen from Fig.QT] only 



7 IKanekar et alj {2007) report a de-convolved source size of 17.6 X 6.6 
mas 2 at 328 MHz. 

8 Many of the source sizes are measured at (VLBI) frequencies, which are 
ofte n many times hig her than that of the redshifted 21-cm line (see table 
2 of lCurranetai]l2005l) . and so in Fig. [T3l (bottom) we have nagged each 
of the DLAs according to how many times larger than 1420/(z a b s + 1) 
the frequency at which the source size is measure d, i.e. tAjize/^bs- The 
recently published low frequency VLBA imaging of lKanekar eta?] (2009a) 
accounts for many of the "green" values, where ^sizc/^bs ^ 3. Since a 
number of these sources appear resolved, we have recalculated the extents , 
rather than using the deconvolved sizes quoted in Kan ekar et al 1 J2009 A . 
which are often smaller than the beam and may be unphysical (e.g. ~ pc 
for 0405-331). 
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Figure 13. Top: The absorber/quasar angular diameter distance ratio versus the absorption redshift. Updated from Curran & Webb 1 2006), where the symbols 
are as per Fig. II II and we include only the DLAs which have been searched to T sp i n /f 100 K (see main text). The iso-redshift curves show how 
DAdla/DAqsq varies with absorption redshift, where DAqsq is for a given QSO redshift, given by the terminating value of z^s (throughout this paper 
we use Ho = 71 km s _1 Mpc -1 , f2 ma ttcr = 0.27 and Q\ = 0.73). That is, we show the DAdla/DAqso curves for z om = 0.5, 1, 2, 3 & 4. Bottom: 
The radio source size versus the absorption redshift. The symbols are colour coded according to the ratio of the frequency of the image to the redshifted 1420 
MHz frequency: green - u s i zc /u^ s 3, orange - 3 < fsize/^abs ^ 10 m ^ re d ~~ "size/fabs > 10. The hatched histogram represents the detections and 
the bold histogram the non-detections. 



three published non-detections are sufficiently sensitive to detect 
r 8pin // > 2000 K at 2 abs > 2 (Sect.O- 

While, due to geometry effects, a high angular diameter dis- 
tance will disadvantage the absorption cross-section in comparison 
to the same absorber placed at a low angular diameter distance, 
with the 21 -cm searches spanning a look-back time range of 12 Gyr 
there are also evolutionary effects to be considered. For instance, an 
evolution in the morphologies of the absorbing galaxies is expected 
to result in different metallicities and these heavy elements can 
provi de cooling pathways for the diffuse gas (Dalgar no & McCravl 
Il972h . Furthermore, it is shown that the relative mix of the cold 
neutral medium (CNM, where T ~ 150 K and re ~ 10 cm -3 ) 
and the warm neutral medium (WNM, T ~ 8000 K and re ~ 0.2 
cm~ 3 ) can be affected by the metallicity (jWol fire et all 1995b . 

At low redshifts, where the DLA hosts can be imaged directly, 
the absorbi ng galaxies appear to be a mix of s pirals, dwarf and 
LSBs (e.g. iLe Bran et al]| 19971 ; iRao et aljr2003l : Ichen & Lanzettd 
120031) and, in general, large galaxies appear to be more popu- 
lous at low redshift with smaller morphologies becoming more 
comm on at high redshifts (Baker et al. 2000l; lLanfranchi & Friacal 
12003b . Therefore, in the more compact absorbing galaxies, where 



the metallicities are lower than in the larger spirals (see Sect. 14.21 ). 
we may e xpect the gas to be m ore dominated by the WNM and, as 
shown bv lCurran et alj ( I2007bl) . T ap i n // is anti-correlated with the 
metallicity (see Sect. l4.2t . However, althou gh the metallicity is it- 
self anti-correlated with the look-back time ( Prochaska et al J2003I : 
IRao et al] |2005h . no [M/HJ-Zaba c orrelation is seen for the 21-cm 
searched sample dCurran et al.l2007bh . and, as seen in Fig.[T3] very 
few morphologies are known for these at z a b s ~ 1. Furthermore, 
the T S pi n /f - [M/H] anti-correlation is consistent with either the 
spin temperature being lower or the covering factor being larger 
at high metallicities, with the larger galaxies likely to have high 
values of / as well as [M/H]. Although the T sp i n // degeneracy is 
unbreakable, it is possible that both factors contribute and that the 
values of [M/H], T sp i n and / are in fact interwoven dCurran et al.l 
l2007bh . 

In addition to decreasing metallicities, the increase in the 
backgro und ultra-violet flux will increase the ionisation fraction of 
the gas (Bajtlik et al. 1988) , in addition to raising its spin temper- 
ature dFieldlll959l ; iBahcall & Ekersi ri969), an effect also seen in 
associated systems, where 21-cm absorption has never yet been de- 
tected in the host galaxy of a quasar with Luv ~ 10 23 W Hz -1 
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Figure 14. In the small angle approximation 8 = c? a bs/ ^^DLA = 
d.QSo/DA QSO , giving d ab s( m in) = 0QSO- DA DLA when the absorber 
just covers the background radio source size. 
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Figure 15. The minimum extent of the absorber required to cover the radio 
source versus the normalised line strength. The symbols are as per Fig. 1131 
(bottom). The line shows the least-squares fit to the detections. The upper 
limits in the ordinate are due to upper limits in the background radio source 
sizes. 



dCurran et alj|2008h . Therefore, as well as always having large an- 
gular diameter distance ratios and lower metallicities, high redshift 
absorption systems will generally be subject to higher ultra-violet 
fluxes, although for intervening absorbers the effect will not be as 
severe as for associated absorbers. Nevertheless, this could be the 
cause of many of the T sp i n // lower limits at z a b s ~ 2.5 (Fig. II lb, 
although, again, such heated gas could have a significantly lower 
21-cm absorption cross-section as well as a higher temperature. 

4.1.2 Absorber extents 

In order to consider the radio source size in conjunction with 
the ratio of the angular diameter distances, we calculate the ex- 
tent of absorber required to fully cover the radio source (Fig. I14t 
and show this against the normalised line strength, 1.823 x 
10 18 / rdv/Nm = //Tspta (Fig. ED. Fr om this we see a corre- 
lation between the extent of the absorber and the 21-cm absorption 
strength for the detections. For the 18 detections this is significant 
at 2.7a, falling to 2.2<r for the 12 detections with J-'size/i-'abs ^ 3 
("green"), with the non-detections exhibiting no correlation what- 
soever (bringing the significance for the whole 39 strong sample 
down to 1.4crO 



Using the deconvolved values o f Kan ekar etal~U 2009a), rather than those 
measured (see footnote[8j, gives 2.5 , 1.9 & 1.3<r, respectively. 



The fit in Fig. Q3] suggests that the non-detections which have 
been searched sufficiently deeply, (i.e. neglecting the two with 
//Tspin ~ 0.01), may require large absorbers (> 10 kpc scales) in 
order to achieve effective covering factors. As is seen from the fig- 
ure, several detections also occupy t his regime, although these ten d 
to be spiral galaxies. Therefore, li ke Kanekar & Cheng alurl d2003l) . 
and as discussed in ICurran et al.l 42005), the majority of non- 
detections (which are at high redshift and ha ve unknown morpholo- 
ies) m ay be non-spirals. Although, unlike IKanekar & Chengalurl 
20031) , we again suggest that may be their small sizes, in addition 
to/rather than high mean spin temperatures, which may be respon- 
sible for their weak line strengths. 

The covering factor can be defi ne d^) as / = 
[d abs /(8Qso-DA DLA )] 2 dCurran & Webbl 120061) . meaning 
that, in general"! the ordinate in Fig.[l5]is equivalent to d a bs/\/7- 
This therefore suggests, at least for the detections, that the line 
strength increases with the size of the 21-cm absorption region, a 
trend which is apparent, even when the absorber size is weighed 
down by the \/J factor, and, in general, from the normalised line 
strength, 



1.823 x 10 1 



Jrdv f 



oc 8qso- -DAdla = 



rf a bs 

7r 



(2) 



That is, d a bs oc f s ^ 2 /T sp m, which makes sense in terms of both 
covering factor and spin temperature - a larger absorber implies 
a larger covering factor, as well as the possibility of a lower spin 
temperature. Presuming that this is not dominated by the numera- 
tor, this could suggest that t he larger absorbing galaxies a re subject 
to lower spin temperatures (Ka nekar & Chengalurll2003l and refer- 
ences therein). However, one interpretation of the distribution of 
the non-detections along the ordinate in Fig. [T5]is that for a given 
absorption cross-section (i.e. assuming the same d a b s for all the 
DLAs), the covering factor is generally lower than for a number of 
the detections. 

Note that, in addition to geometry effects, it is possible that 
these lower covering factors could also be caused by lower beam 
filling factors, introduced by structure ("dumpiness") in the ab- 
sorbing gas, which may be behind the 21-cm variability seen a t 



z a bs = 0.313 towards 1127-145 ( Kanekar & Chengalur 
and at z a bs = 3 . 386 towards 0201+1 13 ideBruynetal 



2001) 



iBriggs etalj|l997l : IKanekar et at] 120071) . iGupta et alj ( 120091) also 



1996 



note from their large survey of Mg II absorption systems that the 
decreasing number density of 21-cm absorbers with redshift runs 
counter to that o f the Lyman-a and W* 2796 > 1 A Mg II absorbers 
jRao et al ] |2006h . This is attributed to a significant evolution in ei- 
ther the CNM filling factor or the covering factor, the former of 
which could be the result of a relatively extensive WNM, thus in- 
dicatin g a higher mean harmonic spin temperature (Kane kar et al.l 
l2009al) . although the larger fractions of warm gas would not ex- 
clude the presence of cold clumps. Either or both of these scenarios 
would be tied up in the T sp i n /f degeneracy. Nevertheless, 21-cm 



10 It can also be estimated as the r atio of the compact unresolved com- 
ponent's flux to the total radio flux (Briggs & Wolfe 1983; Kanekar et al. 
2009a). However, even if high resolution radio images at the redshifted 21- 
cm frequencies are available, this method gives no information on the extent 
of the absorber (or how well it covers the emission). Furthermore, the high 
angular resolution images are of the continuum only and so do not give 
any information about the depth of the line when the extended continuum 
emission is resolved out. 

11 That is, where d abs is the actual extent of the absorber, rather than the 
extent minimum required to cover 8qsq. 
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Figure 16. The normalised line strength (= //T ap i n ) versus the metallic- 
ity for the DLAs searched in 21 -cm absorption, where the symbols are as 
per Fig. 1111 The significance of the correlation is given for both the whole 
sample and the detections only and the least-squares fit for the detections 
shown. 



detection rates are significantly lower at high redshift and the filling 
factor of the CNM, through either its effect on T sp i n or /, could be 
a major contributor to this. 



4.2 21-cm line strength-metallicity correlation 

Addressing the 21-cm line strength- metallicity co rrelation 
(Tspin/. f ~ [M/H] anti-correlation) of ICurran et ail <2007bh . the fo- 
cus of the 2008 GBT observations (Sect. [2~3l l. it appears that we 
should have detected absorption in at least one of the two DLAs for 
which reasonable sensitivities were reached (Fig. |16t . Regarding 
the best improved limi{3, 21-cm absorption may not have been 
detected towards 0528-250, on the grounds that, with an absorp- 
tion redshift very close to that o f the emission redshift (both at 
2 = 2.8110, iMoller et al.l 2004 and references therein), the ab- 
sorber is possibly located in the host of the background quasar. 
We may therefore expect the neutral gas to be subject to excita- 
tion effects seen in other high redshift 21-cm absorption searches 
dCurran et al1l2008h . However, the presence of H2 in the absorber 
suggests that the gas is relatively cool with T^ n w 110 — 150 K 
(Srianand et al. 2005) and so this is unlikely. 

In any case, referring to Fig. [T6j the limit is not overly low 
in comparison to 0438-436 (the detection at [M/H] = -0.68, 
//Tspin = 6.3 x 10~ 4 [KT 1 ]) and there is also the possibility that 
we have missed the line, with our limit only being good over the 
redshift range 2 w 2.806 - 2.812 (Sect.|23l). Assuming this limit 
is reliable and inc luding 0311+430, strengthens the correlation (cf. 
ICurran et all2007bh. 

lYork et al.l fe007h comment that the 2 a b s = 2.289 absorber 
towards 0311+430 has an unusually high metallicity, although ac- 
cording to Fig. [16] this has some leeway before the metallicity 
becomes atypically high. From its grouping in the line strength- 
metallicity pl ot, we suggest th at the absorption is due to a spiral 
galaxy, as do I York et al. J2007h on the basis of its low "spin tem- 
perature". Support for this assertion is given by its grouping with 



the spirals in Fig. 1151 (the "unknown" at 7 x 10 3 K 1 , 5 kpc), 



12 By a factor of « 7, cf. lCarilli et aljl99fj . See footnote^] 



although this is on the grounds of a large absorption cross-section, 
rather than the spin temperature. 



5 SUMMARY 

Through three sets of separate observations, we have undertaken 
a survey for 21-cm absorption in a total of twelve DLAs and sub- 
DLAs. Due to RFI and low neutral hydrogen column densities, use- 
ful sensitivities were obtained for five of the sample, none of which 
resulted in a detection. We add these to the recent results from other 
surveys, which while forming a very heterogeneous sample, sug- 
gest that, at least statistically: 

(i) There appears to be no large difference in the spin tempera- 
ture/covering factor ratio between the low and high redshift DLAs 
(Fig.EUand Tabled, although the non-detections do have a higher 
mean redshift and more sensitive observations of these would be 
required to fully address this. Nevertheless, the fact that the mean 
Tspin// difference is largest for a redshift partitio n of z^ -,^ = 1.0 



is con sistent with the geometry effects noted by ICurran & Webbl 
(2006). Other factors, such as increasing WNM fractions (increas- 
ing Tspin) or an evolution in morphologies (also decreasing /) with 
redshift cannot be ruled out, although the T sp i n // degeneracy can- 
not be broken. 

(ii) The observations since ICurran et al] J2005h confirm that the 
21-cm line strength is correlated with the total neutral hydrogen 
column density, which further suggests that the spin temperature 
to covering factor ratio (T sp i n //) is roughly constant, for both the 
detections and non-detections (Fig. 1121). 

(iii) By combining the radio source sizes with the angular di- 
ameter distances to the DLAs, we show that the absorber extent 
normalised by the covering factor (d a bs/\/7) is correlated with the 
line-strength for the detections, with the non-detections generally 
exhibiting larger values of d„hs/\/T (Fig.ll5t. A possible interpre- 
tation is that the non-detections require larger minimum absorber 
extents in order to be detected. This may also be interpreted as the 
covering factors being lower than for the detections. 

This argument relies on the assumption that the unknown 
morphologies of the non-detections are not those of spi- 
ral galaxies. This is consi stent with both t he hy potheses of 
Kane kar & Chen galur (200^) and lCurran etai] [2005). in that spi- 
ral galaxies may have lower spin temperatures or larger covering 
factors than the other morphologies, respectively. The constancy of 
Tspin// between the detections and non-detections (point ii, above) 
may suggest that both contribute, with lower covering factors be- 
ing balanced by lower spin temperatures. The absorber extent-line 
strength correlation further suggests that d a b s oc / 3 ^ 2 /T sp in for the 
detections. Although, once again, we are limited by the spin tem- 
perature/covering factor degeneracy, we do believe that the large 
scatter in this (T sp i n / / = 60 to > 13 000 K) is mainly dominated, 
or at least balanced, by different covering factors: This represents 
a vast range of spin t emperatures (specifically, T sp in ~ 200 to > 
9 000 K suggested bv lKanekar & Chengalurll2003h . Although hav- 
ing a much smaller dynamic range than the spin temperature, the 
covering factor is subject to three separate effects - the extent of 
the absorbing region, the size and morphology of the background 
1.4 GHz emitter, as well as the effects introduced by the geome- 
try of a flat expanding Universe. Furthermore, the beam filling of 
the cold neutral medium in these absorbers could have a significant 
effect on both the covering factor and th e spin temperature . 

(iv) We agree with the suspicions of lYork et alj d2007l) that the 
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z a bs = 2.289 absorber towards 0311+430 arises within a spi- 
ral galaxy. Although not on the basis of its low spin temperature 
(which we do not find unusual at these redshifts), but on its ab- 
sorption cross-section and its grouping with the other spirals on a 
line-strength-metallicity plot. 
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